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Abstract. Let fi C R , d ^ 1, be a bounded domain with piecewise 
smooth boundary dft and let U be an open subset of a Banach space 
Y . Motivated by questions in "Uncertainty Quantification," we consider 
a parametric family P = (P y ) y eu of uniformly strongly elliptic, second 
order partial differential operators P y on Q. We allow jump discontinu- 
ities in the coefficients. We establish a regularity result for the solution 
u : O x U — > K of the parametric, elliptic boundary value/transmission 
problem P y u y = f y , y € U, with mixed Dirichlet-Neumann bound- 
ary conditions in the case when the boundary and the interface are 
smooth and in the general case for d = 2. Our regularity and well- 
posedness results are formulated in a scale of broken weighted Sobolev 
spaces 1C™+l{Q) of Babuska-Kondrat'ev type in Q, possibly augmented 
by some locally constant functions. This implies that the parametric, 
elliptic PDEs (P y ) y su admit a shift theorem that is uniform in the 
parameter y £ U. In turn, this then leads to /i m -quasi-optimal rates 
of convergence (i. e., algebraic orders of convergence) for the Galerkin 
approximations of the solution u, where the approximation spaces are 
defined using the "polynomial chaos expansion" of u with respect to a 
suitable family of tensorized Lagrange polynomials, following the method 
developed by Cohen, Devore, and Schwab (2010). 

1 Introduction 

Recently, questions related to differential equations with random coefficients have 
received a lot of attention due to the practical applications of these problems 
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[2 7 9116117] . Our paper is motivated by the approach in [6113] , where families 
of differential operators on polyhedral domains indexed by y € U, were studied. 
As in those papers, U is an open subset of a Banach space Y , which allows us 
to study the analyticity of the solution in terms of y G U . 

We here study the properties of solutions to a family of strongly elliptic, 
mixed boundary value/transmission problems 

PyU y (x) = Pu(x,y) = f y (x) = f(x,y), x€f2,y€U (1) 

on a domain Q C M. d , d > 1. The domain SI is assumed to be piecewise 
smooth and bounded. Thus, for each y g U, we are given a second order, 
uniformly strongly positive, parametric partial differential operator P y on SI 
whose coefficients are functions of (x, y) 6 Q x U and are allowed to have 
jump discontinuities across a fixed interface r. More precisely, we assume that 
Q = U^ =1 f2k, where Slk are disjoint domains with piecewise smooth boundaries 
and r := ( uf =1 dQ k ) \ dSl. 

Under suitable regularity assumptions on the coefficients of P and on the 
source term / : fi x U — > K, we establish in Section [5] a regularity and well- 
posedness result for the solution u : SI x U —> M of the parametric, elliptic bound- 
ary value/transmission problem ((TJ with mixed Dirichlet-Neumann boundary 
conditions. Our regularity result is formulated in a scale of broken weighted 
Sobolev spaces t^(Q) = ®J^ 1 K£fi L (Sl k ) of Babuska-Kondrat'ev type in SI, 
for which we prove that our elliptic PDEs (P y ) ye u admit a shift theorem that 
is uniform in the parameter y S U. We deal completely in this paper with the 
cases when the boundary dS7 and the interface -T are smooth and disjoint. We 
also indicate how to proceed in the general case for d = 2. Our results gener- 
alize the results of |13j by allowing jump discontinuities in the coefficients and 
by allowing adjacent edges to be endowed with Neumann-Neumann boundary 
conditions. We will be therefore brief in our presentation, referring to [13] . as 
well as [61 7) for more details. 

The main contribution of this paper is to study the regularity of the solution 
of a (non-parametric) transmission/boundary value problem with rather weak 
smoothness assumptions on the coefficients. As far as we know, this paper is the 
only place where a complete proof for the regularity of transmission problems is 
given, even in the case of smooth coefficients. The results are general enough so 
that one can use the approach in |6ll3j to obtain regularity results for families 
and then to obtain optimal rates of convergence for the Galerkin method. An 
abstract version of this method is explained in [3] . These issues will be discussed 
in more detail in a forthcoming paper. 

The paper is organized as follows. In Section [5] we formulate our parametric 
partial differential boundary value/transmission problem and introduce some 
of our main assumptions. We also discuss the needed notions of positivity for 
families of operators and derive some simple consequences. In Section [3l we 
review the "broken" version of usual Sobolev spaces, and then formulate and 
prove the main results, Theorem [TJ which is a regularity and well-posedness 
result for non-parametric solutions in smooth case. In Section 0] we recapitulate 
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regularity and well-posedness results for the non-parametric, elliptic problem 
from [415110] . the main result being Theorem[2J This theorem is then generalized 
to families in Section [5j thus yielding our main regularity and well-posendess 
result for parametric families of uniformly strongly boundary value/transmission 
problems, namely Theorem [31 As mentioned above, this result is formulated in 
broken weighted Sobolev spaces (the so called "Babuska-Kondrat'ev" spaces). 
As in |6I7I13) these results lead to /i m -quasi-optimal rates of convergence for a 
suitable Galerkin method for the approximation of our parametric solution u. 

Acknowledgements. VN acknowledges the support of the Hausdorff Institute 
for Mathematics (HIM) in Bonn during the HIM Trimester "High dimensional 
approximation," where this work was initiated. We also thank Christoph Schwab 
for several useful discussions. 

2 Ellipticity, positivity, solvability for parametric families 

We now formulate our parametric partial differential boundary value/transmission 
problem and introduce some of our main assumptions. 

2.1 Notation and assumptions 

By Q C M d , d > 1, we shall denote a connected, bounded piecewise smooth 
domain, which we assume is decomposed into finitely many subdomains flk 
with piecewise smooth boundary, (7 — L)^ =1 f2k- We obtain results on the spatial 
regularity of PDEs whose data depend on a parameter vector y £ U C Y , where 
U is an open subset of a Banach space Y. By a l J q ,b pq ,c pq : fl x U — >• M, 1 < 
hi _• d, we shall denote bounded, measurable functions satisfying smoothness 
and other assumptions to be made precise later. We denote by A = (a^ q , b pq , c pq ). 
Let us denote by 9, = i — 1, . . . , d. We shall then denote by P A — [P pq ] a 
fi x fi matrix of parametric differential operators in divergence form 

(d d \ 

-^2d l (a p ] q (x,y)d j ) + '*Tb t pq (x, y)d l + c pq (x, y) u(x,y), (2) 
i,j=i i=i / 

where x £ fl and y E U . Note that the derivatives act only in the x-direction 
and y is just a parameter. The matrix case is needed in order to handle the case 
of systems, such as that of (anisotropic) linear elasticity. 

A matrix P = [P pq ]p ?=1 of differential operators acts on vector-valued func- 
tions u = {u q ) q=1 in the usual way (Pu) p — J2q=i Ppq u qi f° r u = { u q) 
C°°(n x uy. We recall that H- 1 {Q) is defined as the dual of H^{Q) := {u e 
H 1 (S7), u\oa = 0} with pivot L 2 (f2). Occasionally, we shall need to specialize a 
family P for a particular value of y, in which case we shall write P y : C°°(f2) fl —> 
i? _1 (J7) M for the induced operator. We emphasize that we allow P to have non- 
smooth coefficients, so that Pu may be non-smooth in general. 
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2.2 Boundary and interface conditions 

We impose mixed Dirichlet and Neumann boundary conditions. To this end, we 
assume there is given a closed set do^ C df2, which is a union of polygonal 
subsets of the boundary and we let dj^fl := dfl s do^- The set do^ will be 
referred to as "Dirichlet boundary" and d/vf? as "Neumann boundary," accord- 
ing to the type of boundary conditions that we associate to these parts of the 
boundary. The case of cracks is also allowed, provided that one treats different 
sides of the crack as different parts of the boundary, as in 10 , for instance, 
but we choose not to treat this case explicitly in this paper. We then define the 
conormal derivatives 

ju d 

(V A u) p (x, V) = ^2 ^2 Via.p J q (x, y)dju q (x, y), x G d N Q, y G U, (3) 

where v = (yf) is the outward unit normal vector at x 6 dtqfl. The conormal 
derivatives V^u^ at the interface F are defined similarly, using an arbitrary but 
fixed labeling of the two sides of the interface into a positive and a negative part. 

We shall also need the spaces Hjj(Q) and H^^Q) for vector- valued func- 
tions: H}>(f}) := {u S H x (fiY, u = on d D f2} and H] D 1 {fl), defined to be the 
dual of Hjj(f2) with pivot space L 2 (fl). Note that we assume here as in [13] 
that we have the same type of boundary conditions for all components u q of the 
solution vector u. 

Recall that our domain fl is decomposed into K subdomains of the same 
type (with piecewise smooth boundary), fl = U^ =1 f2k- We then denote by r := 
( ^k=i x ®^ tli e interface of our problem. We also fix arbitrarily the sides 
of the interface, and we thus denote by u + , respectively u~ the non-tangential 
limits of u at the two sides of the interface. We define similarly the conormal 
derivatives V^u + and V^w~ at the interface, but using the two sided limits 
of the coefficients a^ q at r. We consider the parametric family of boundary 
value/interface problems 

P A u(x,y) = f(x,y) x G Q, 

u(x,y) = x £ dofl, 

V^u(x, y) = g{x, y) x € d N f2 (4) 

u + (x, y) — u~(x, y) = ief 

V A u + (x, y) — V A u~ (x, y) = h(x, y) ief 

where P A is as in Equation ^j, \7 A is as in Equation (O, and y G U. We stress 
that for us the dependence of P A on its coefficients, that is on A, is important, 
which justifies our notation. 



2.3 Ellipticity and positivity for differential operators 



In this subsection we recall the definition of the positivity property for parametric 
families of differential operators. Let us therefore consider, for any y G U, the 
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parametric bilinear form B(y; • , • ) defined by 

B(y;v,w) := J xen T, P , q =i {T,i,j=i a '%(x,y)diV p (x,y)d j w q (x,y) 

T,i=ib pq (x,y)diV p {x,y)w q (x,y) + c pq (x,y)v p (x,y)w q (x,y)jdx , y EU. (5) 

Definition 1. The family (P y ) y ^u is called uniformly strictly positive definite 
on Hq ( Q)^ cl'C H 1 if the coefficients a l fj q are symmetric ini,j and in p, q 
(that is, a p J q = a p l q = a q J p , for all i, j, p, and q), and if there exist < r < R < oo 
such that for all y E U , and v, w E V , we have 

\B(y;v,w)\ < R\\v\\ H i (n) \\w\\ H i {n ) and r\\v\\ 2 H1[n) < B{y;v,v). 

If U is reduced to a single point, that is, if we deal with the case of a single 
operator instead of a family, then we say that P is strictly positive definite. 
Throughout this paper, we shall assume that {P y ) y ^u is uniformly strictly pos- 
itive definite. Positivity is closely related to ellipticity. 

Definition 2. The family {P y ) ye .jj is called uniformly strongly elliptic if the 

coefficients a p 3 q are symmetric in i,j and in p, q and if there exist < r e < R e < 
oo such that for all x E D, y EU , £, E R d , and rj E W 

ft d 

r e \e\rf < E E <&(*>v)ZiZjW 9 < Re\e\v\ 2 ■ (6) 
p,g=i »,i=i 

In case one is interested only in scalar equations (not in systems), then for 
V = Hj D (Q), the assumption that our family P y is uniformly positive definite 
can be replaced with the (slightly weaker) assumption that the family P y is 

uniformly strongly elliptic, that Yli=i — m Ei=i — on c^i?, 
c > 0, and dr>Q ^ (in which case it also follows that SdJ? has a non-empty 
measure). In general, a uniformly strictly positive family P will also be uniformly 
strongly elliptic. 

2.4 Consequences of positivity 

The usual Lax-Milgram lemma gives the following result as in [6]. Recall the 
constant r from Definition [1] 

Proposition 1. Assume that f y := f(-,y) E H^ ) 1 {Q), for any y E U. Also, 
assume that the family P is uniformly strictly positive definite. Then our family 
of boundary value problems P y u y = f y , u y E Hp(fi), i.e., Equation O), admits 
a unique solution u y — P y x f y . Moreover, \P y \h{h^-h\ ) — r > f or a ^ V e 

The parametric solution u y E Hj-,(0) of Proposition Q] is then obtained from 
the usual weak formulation: given y E U , find u y E V := Hp(f2) such that 

B(y;u y ,w) — (f y ,w) + / g y wdS + / h y wdS, Vw E V, (7) 
Jd N n Jr 

where {f y ,w) denotes the L 2 {f2) inner product and dS is the surface measure 
on dfl or on P. Also, f y (x) = f{x,y), and similarly for u y , g y , and h y . 
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3 Broken Sobolev spaces and higher regularity of 
non-parametric solutions in the smooth case 

One of our main goals is to obtain regularity of the solution u both in the space 
variable x and in the parameter y. It is convenient to split this problem into two 
parts: regularity in x and regularity in y. We first address regularity in x in the 
case when the boundary dfl and the interface P are smooth and disjoint. We also 
assume that each connected component of the boundary is given a single type 
of boundary conditions: either Dirichlet or Neumann. This leads to Theorem [TJ 
which states the regularity and well-posedness of Problem ^ in this smooth 
case (dfi and F smooth and disjoint). This is the main result of this paper, 
and, as far as we know, no complete proof was given before. We also consider 
coefficients with lower regularity than it is usually assumed, which is needed to 
treat the truly parametric case. (We are planning to address this question in a 
future paper.) 

We assume throughout this and the following section that we are dealing with 
a single, non-parametric equation (not with a family), that is, that U is reduced 
to a single point in this subsection. We also assume that ji = 1, to simplify the 
notation. 

We shall need the "broken" version of the usual Sobolev spaces to deal with 
our interface problem. Recall the subdomains fl k C f2, 1 < k < K, we define: 

H m {Q) := {v : Q -> R,v G H m {Q k ), VI < k < K} (8) 
W m '°°{Q) := {v : Q -> M, d a v G L°°(Q k ), \/\<k<K, \a\ < m}. 

For further reference we note that the definitions of these spaces imply that 
the multiplication and differentiation maps W m,oa {Q) x H m (Q) — > H m (il) and 
d % : H m {f2) -4 H m - 1 {f2) are continuous. 

One of the difficulties of dealing with interface problems is the more compli- 
cated structure of the domains and ranges of our operators. When m = 0, we 
define V m = V = H^(f2) = V and K m = TZ = H^(O) = V*. Then we define 
P A in a weak sense using the bilinear form B introduced in Equation ([5]) (see 
the discussion around Equation (2.12) in [TU] for more details or the discussion 
around Equation (20) in [H]). Assume now that m > 1. We then define 

V rn := H m+1 (Q) n {u = on d D Q} n {u+ - u~ = on r} and 

K m := e ff m - 1/2 (9jv^) © H m - l ' 2 {r). 

In particular, V m = H rn+1 {Q) n H^Q). Let A = {a 11 ,b\ c) G W m >°° 

and P A u = Y^i j=i di{a lJ dju) + Y^i=i ° l diU + cu, as before. Then the family of 

partial differential operators P A : V m — } 1Z m , 

P„ A u = (P A u,V A u\ dNn , (V A u+ - V A u+)\ r ) (9) 

is well defined. Note that the domain T> m and codomain lZ m are independent of 
y G U, which justifies why we do not consider homogeneous Neumann boundary 
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conditions. We are now ready to state and prove our main theorem. Let us denote 

( k \ ^ 2 K 

IMIh>"(^) := [T, k =i \\u\\H^(n k y) and |M]^ m ,«> (fi) := E fe =i ll»||w"».«>(fl fc ) 

the resulting natural norms on the spaces introduced in Equation ([5]). 

Theorem 1. Let us assume that £1 C R d is smooth and bounded, that the inter- 
face r is smooth and does not intersect the boundary, and that to each compo- 
nent of the boundary it is associated a single type of boundary conditions ( either 
Dirichlet or Neumann). Assume that A = (a y , c) G W m '°°{Q) d2+d+1 and that 
P A is strictly positive definite on Hjj(il), then P A : T> m — > 1Z m is invertible. 

Moreover, let \\P^ 1 \\ denote norm of the inverse of the map P : Hp(f2) — > 
Hjj{Q)* =: H^ D 1 (n). Then there exists a constant C\ > such that the solution 
u of satisfies 

\HH^ { n) + \Hm { n) < C'i(ll/ll^-i(«)ll5ll Jfm -i (ajv „ ) + ll^ll ifm -i (r) ). (10) 
with the constant C\ — Ci(to, ||P _1 ||, ||A||^ /m =0 ). 

Proof. In the case of the pure Dirichlet boundary conditions for an equation and 
without the explicit bounds in Equation (|10j) . this lemma is a classical result, 
which is proved using divided differences and the so called "Nirenberg's trick" 
(see [8112) ) . Since we consider transmission problems and want the more explicit 
bounds in the above Equation (ITU1) . let us now indicate the main steps to treat 
the interface regularity following the classical proof and [12]. The boundary 
conditions (1 e., regularity at the boundary) were dealt with in |13] . In all the 
calculations below, all the constants C in this proof will be generic constants 
that will depend only on the variables on which C\ depends (i. e., on the order 
m, the norms ||-P _1 || and mil^m.oo^))- We split the proof into several steps. 

Step 1. We first use Proposition Q] to conclude that P : Hj->(f2) — > H]j(Q)* is 
indeed invertible. This provides the needed estimate for to = (in which case, 
we recall, our problem Q has to be interpreted in a weak sense). 

Step 2. For to > we can assume g = and h = by using the extension 
theorem as in [T5] . 

Step 3. We also notice that, in view of the invertibility of P for to = and 
since u 6 H^Q), it suffices to prove 

\\uU m+1[n) < c(J2 II/IIh-^) + IMIff-W • (ii) 

k 

Indeed, the desired inequality (|10p will follow from Equation (jlip by induction 
on m. Since Equation (fTTj) holds for P if, and only if, it holds for A + P, in order 
to prove Equation (fTT|) , it is also enough to assume that A + P is strictly positive 
for some A G M. In particular, Equation (fTTj) will continue to hold-with possibly 
different constants-if we change the lower order terms of P. 

Step 4- Let us assume that f2 — M. d with the interface given by r = {xd = 0}. 
Let fl\ = Ml and fli = be the two halves into which K d is divided (so 
K = 2). Then we prove Equation for these particular domains and for 
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g = and h — by induction on to. As we have noticed, the Equation is 
true for to = 0, since the stronger relation (|10D is true in this case. Thus, we 
shall assume that Equation (1111) has been proved for to and for smaller values 
and we will prove it for to + 1. That is, we want to prove 

||w||^ m +2( R d) < C(||/||^- m(Rd) + ||u||^- m+1(Rd) ). (12) 

To this end, let us first write 

d 

|M|#m + 2( R d) < ^2 ll 5 J U H_ff'« + l(K ti ) + H U lli 2 (K d ) ■ ( 13 ) 

i=i 

We then use our estimate (fTTj) for to (using the induction hypothesis) applied 
to the function djU for j < d. This gives 

II^J U lld m + 1 (R< i ) - \\P^i u \\H^-^(S_d) < WdjfWjjm-i^ + ||[-P)9j]u||^ m _i( R<i j 

<\\f\\ Amm + C\\u\\ 6m+im (14) 

since the commutator [P, dj] = Pdj — djP is an operator of order < 2 whose 
coefficients can be bounded in terms of ||A||^ mi00( -^. We now only need to 
estimate , we do that on each half subspace. 

\\d d u\\fj m+1{ud) < J2 d 3 =i \\9jd d u\\ 6m{m + \\d d u\\ L 2 {md) 

The right hand side of the above equation contains only terms that have already 
been estimated in the desired way, except for ||c^u||#m. Since to > 0, we can 
use the relation Pu = f to estimate this term as follows. Let us write Pu = 
di(a 13 dju) + J2 bl diU + cu. This gives a dd dju = f -J2{i,j)^(d,d) a lJ did-jU + Qu, 
where Q is a first order differential operator. Next we notice that a dd is uniformly 
bounded from below by the uniform strong positivity property (which implies 
uniform strong ellipticity): (a dd )~ 1 < r _1 . Note that by Proposition [TJ we have 
ll-P -1 II < r i an d hence r" 1 is an admissible constant. This gives d\u — (a dd ) -1 /— 
J2j=i B J djU + Q\u where B° and Qi are first order differential operators with 
coefficients bounded by admissible constants, which then gives 

d-i 

\\dd u \\H^(M. d ) - C(ll/ll_H-™(E' i ) H^i M lli?'"+ 1 (K ti ) + ll U ll#"*+ 1 (R d )) ( 16 ) 

J'=l 

< C(ll/llfl-"'(K £i ) + H U Hif m + 1 (R ti )) 

by Equation (|1"3|). Equation (fI5]l and (fTB)) then give 

\\d d u\\ Am+X{VLi) < C(||/|| dm(R<i) + ||u|| dm+ i (R(J) ) . (17) 

Combining Equations (TT7)) and (TT4")) with Equation (fT5|) gives then the desired 
Equation (TT2l for to replaced with m + 1. 
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Step 5. We finally reduce to the case of a half-space or a full space using 
a partition of unity as in the classical follows. We choose a smooth 

partition of unity on SI consisting of functions with small supports. The 
supports should be small enough so that if the support of 4>j intersects the 
boundary of SI or the interface r, then the boundary or the interface can be 
straightened in a small neighborhood of the support of <pj . We arrange that the 
resulting operators are positive and we complete the proof as in |13j . 

See also [THIS] . 

4 Weighted Sobolev spaces and higher regularity of 
non-parametric solutions 

We now assume d = 2, so SI is a plane domain. We allow however f? to be 
piecewise smooth. We also consider coefficients with lower regularity than the 
ones considered in (TTJ] . This leads to Theorem [51 which will be then generalized 
to families in a forthcoming paper, which will contain also full details for the 
remaining results. We continue to assume that we are dealing with a single, 
non-parametric eguation and that \i = 1. 

To formulate further assumptions on our problem and to state our results, 
we shall need weighted Sobolev spaces, both of L 2 and of L°° type. Let V be 
the set of singular points, where Q € V if one of the following is satisfied: (1) it 
is a vertex, (2) it is a point where the boundary condition changes type (from 
Dirichlet to Neumann), (3) it is a point where the interface meets the boundary, 
or (4) it is a non-smooth point on the interface of the subdomains Qk- Let us 
denote from now on by p : R 2 — > [0, 1] a continuous function that is smooth 
outside the set V and is such that p(x) is equal to the distance from x € R 2 to V 
when x is close to the singular set V. The function p will be called the smoothed 
distance to the set of singular points. We can also assume || Vp|| < 1, which will 
be convenient in later estimates, since it will reduce the number of constants (or 
parameters) in our estimates. We first define the Babuska-Kondrat'ev spaces 

:= {v : Q -> R, p^ a d a v G L 2 (f2), V \a\ < m} (18) 
W m '°°(J?) := {v : Q -> R, p^d a v e L°°(f2), V \a\ < to}. (19) 

We shall denote by || • ||/c m (j?) and || ■ ||w" oo (j7) the resulting natural norms 
on these spaces. We shall need also the "broken" version of these Babuska- 
Kondrat'ev spaces for our interface problem. Recall the subdomains ilk C 
1 < k < K . In analogy with the smooth case, we then define: K™{ST) := {v : 
Q -> R, v e fC™(n k ), VI < k < K}, and VV m ^°°(^) := {v : Q -> R, v G 
W" h00 (f2 k ), VI < k < K}. If V is empty (that is, if the domain Q is smooth 
and the interface is also smooth and does not touch the boundary), then we 
set p = 1 and our spaces reduce to the broken Sobolev spaces H m (f2) and 
W m ^°° introduced in the previous section, Equation ([S]). As in the smooth case, 
the multiplication and differentiation maps VV m -°°(J?) x K™(fl) -> K™{ti) and 
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di : K.™{Q) — > K™~i{fi) are continuous. Let S C dfik- Also as in the smooth 
case, we define the spaces '{S) as the restrictions to S of the functions 

u e ]C™+i(f2k)- These spaces have intrinsic descriptions [llllj similar to the 
usual Babuska-Kondrat'ev spaces. Note that no "hat" is needed for the boundary 
version of the spaces t. Also K%£fi(Si U S 2 ) = (Si) © K^ff(S 2 ), if 

Si and S 2 are disjoint. 

We need to consider the subset V s of V consisting of Neumann-Neumann 
corners (1 e., corners where two edges endowed with Neumann boundary con- 
ditions meet) and non-smooth points of the interface, which can be described 
as V s := V n {Q £ V, Q 6 doL?}. Note that, if a point Q at the intersection 
of the interface F and the boundary falls on an edge with Neumann boundary 
conditions, then Q is also included in V s . In order to deal with the singularities 
arising at the points in V s (which behave differently than the singularities at the 
other points of V), we also need to augment our weighted Sobolev spaces with 
a suitable finite-dimensional space. Namely, for each point Q £ V s , we choose a 
function xq S C°°(f2) that is constant equal to 1 in a neighborhood of Q. We 
can choose these functions to have disjoint supports. Let W s be the linear span 
of the functions XQ f° r anv Q £ V s . We now define the domains and ranges of 
our operators. Assume first that m > 1. 

V a>m := (K%fi(£!) + W s ) n {u = on d D f2} n {u+ - u~ = on T} 

n a , m ■.= t^{Q) © K™:t/*{d N n) © K-i/'in- 

Let us observe that, by definition, the functions in W s satisfy the interface and 
boundary conditions (so W s C V := Hj } (f2)). Moreover, for a > 0, we have 
D om = (K^i(n) + W a )nHjj(n). Denote A = (a y ',6 l ,c) G w™^(f2) d2+d +\ 
d — 2, and P A u = Y^V . =1 d^a^djU) + Yli=ib z diU + cu, as before. Then the 
family of partial differential operators P^ m : X> a , m — > 7tL a]m 

P a A mM = (Pu^u\ dN n, (V;V - V^+)| r ) (20) 

is well defined and the induced map w™,°°(Q)(d 2 +d+i) 3 A = (a ij ,b\c) -> 
S £(2? , m , T^ a ,m) is continuous (recall that d = 2). The continuity of this 
map motivates the use of the spaces VV m '°°(J7). 
When TTi — 0, we define 

Da,m = ^a,o = Ka+i( fi ) H {m = on d D .f2} + W s 
K a , m = n afl = (K}_ a+1 (Q) n {u = on dptf})*, 

where in the last equation the dual is defined as the dual with pivot L 2 (fl). Then 
we define P aj o in a weak sense using the bilinear form B introduced in Equation 
([5]), as in the smooth case. 

Recall the constant < r in the definition of the uniform strict positivity 
(Definition [T]) and Proposition [TJ We now state the main result of this section. 
Recall that U is reduced to a point in this section. 
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Theorem 2. Assume that A = (a^,b l ,c) G yV m '°°(r2) <i2 + << + 1 J d = 2, and that 
P A is strictly positive definite on V = (.!?). Then there exists < r\ such 
that for any m G No and for any < a < rj, the map P A m : T> a m — > TZ a , m is 

boundedly invertible and ||(-fX m ) -1 || — where C — C(m, r, a, ||^-Hvv"».°°(j7)) 
depends only on the indicated variables. 

A more typical formulation is given in the following corollary. 
Corollary 1. We use the notation and the assumptions of Theorem^ If f G 
9 e K-™-i/2 (9n^), andh G K^lj/jf (-Oj then the solution u G H})(f2) 
of Problem can be written u = u r + u s , with u r G K7£^f(f2) and u G W s , 
such that 

\\u r \\ t7 +i + \\u s \\ L 2 < (7(11/11^-1 + llsllx^_-y/( Sjv i3) + H^Hx^-y/cr))' 
with C as in Theorem^ 

5 Applications 

We keep the settings and notations of the previous section. In particular, d = 2 
and we are dealing with equations (not systems) . One can proceed as in [617113] to 
obtain /i m -quasi-optimal rates of convergence for the Galerkin u n approximations 
of u. Namely, under suitable additional regularity in the y G U variable one can 
construct a sequence of finite dimensional subspaces S n C L 2 (U; V) such that 

\\u-u n \\ L 2 (u . v) <Cdim(S n )- m / 2 \\f\\ Hm - Hn) . (21) 

This is based on a holomorphic regularity in U and on the approximation prop- 
erties in [10] . We now state a uniform shift theorem for our families of boundary 
value/transmission problems. 

Let us denote by C£(U; Z) the space of fc-times boundedly differentiable func- 
tions defined on U with values in the Banach space Z. By C£(U;Z) we shall 
denote the space of analytic functions with bounded derivatives defined on U 
with values in the Banach space Z. Recall that r is the constant appearing in 
the definition of uniform positivity of the family (P B ) Bg [j. Theorem [5] extends to 
families of boundary value problems as in |13| as follows. Let us denote by rj(y) 
the best constant appearing in Theorem [2] for P = P y and r/ = mi ye u r](y). 

Theorem 3. Let m G No and fco G No U {oo, ui} be fixed. Assume that A = 
{a %1 ,b\c) G C^ (U;W m ' oa (Q)) d2+d+1 ; d = 2, and that the family P A is uni- 
formly positive definite. Then r\ — inf ye jj rj(y) > 0. Let f G C b ° {U^KJ^Zx {fi)), 
g G Cl {U;K^Ii/2{d N ^)), h G tf°(U; K^~^{r)), and < a < r). Then the 
solution u of our family of boundary value problems ^ satisfies u G C b ° (U; £> a ,m)) 
Moreover, for each finite k < ko, there exists a constant C a , m > such that 

\\ U \\c*{U;V a , m ) < ^a,m(||/|| C k(y. 

+ \\9\\ c k {u . K ™-ij2 {dN[1)) + \\h\\ c k {u . K ™-i/2 {r)) ) . 
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Hengguang Li, Victor Nistor, and Yu Qiao 



The constant C ajTn depends only on r, m, a, k, and the norms of the coefficients 
a ij ,b\c in C^{U; W m >°°(f2)), but not on f or g. 
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